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Abstract 

We study the effects of threshold resummation on the rapidity dependence of single- 
particle-inclusive cross sections, using the prompt photon cross section as an example. We 
make use of the full resummation formula at next-to-leading logarithmic accuracy and de- 
velop a new technique for treating rapidity in resummation. We compare our phenomeno- 
logical results with those of previous studies and discuss differences and similarities of the 
two existing resummation formalisms. 
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1 Introduction 



Prompt-photon production at high transverse momentum [|IJ, pp,pp,pN — > 7X, has been a clas- 
sic tool for constraining the nucleon's gluon density, because at leading order a photon can be 
produced in the Compton reaction qg — > ^yq. The "point-like" coupling of the photon to the 
quark provides a potentially clean electromagnetic probe of QCD hard scattering. A pattern 
of disagreement between theoretical predictions and experimental data for prompt photon pro- 
duction has been observed in recent years |2|, [3[ |], however, not globally curable by changing 
the factorization and renormalization scales in the calculation, or by "fine-tuning" the gluon 
density §, 0- The most serious problems relate to fixed-target data, where next-to-leading 
order (NLO) theory dramatically underpredicts some data sets f|. At collider energies M, 
there is less reason for concern, but also here the agreement is not fully satisfactory. In this 
context, threshold resummations for the inclusive prompt photon cross section have been devel- 
oped in [^, ^ and applied phenomenologically in [ ffi, |T1[] . At partonic threshold when the initial 
partons have just enough energy to produce the high-p T photon and a massless recoiling jet, the 
phase space available for gluon bremsstrahlung vanishes, resulting in logarithmic corrections to 
the partonic cross section. Threshold resummation || ^ [H| |T^| organizes these corrections to 
all orders in a,. 



Phenomenological studies |TQ, 11] of threshold resummation for the prompt photon cross 



section have shown that the size of the resummation effects increases with xt = 2pt/VS, where 
Pt is the photon transverse momentum. This is expected || [| from the interplay of the resummed 
partonic cross section with steeply falling parton distributions, and explains why resummation 
effects are small for the prompt photon cross section measured at the Tevatron |2[], where xt is 
generally smaller than 0.1. In contrast to this, for fixed-target energies, where values as large 
as xt ~ 0.7 are attained 0, a significant enhancement of the theory prediction was found. 
The enhancements found in |[10| , |TT|j , however, are not sufficient to bring theory in line with all 
fixed-target data (or, of course, to reconcile conflicting data sets). Nonetheless, valuable insights 
have been gained from applying threshold resummation to the prompt photon cross section. 
One of them is that after resummation the theory predictions show a significant improvement 
in factorization and renormalization scale dependence |T(| over NLO, which has a notoriously 
large scale sensitivity. This implies that threshold resummation does provide an important 
addition to the theory result, and has helped develop approaches |L4|, [L5| that extend threshold 
resummation, by redistributing its contributions to the hard scattering function over soft gluon 
transverse momenta, giving rise to further enhancement of the cross section through recoil effects 
on the prompt photon spectrum. 

We believe that further analyzing the effects of pure threshold resummation is useful, for 
two reasons. First, as mentioned above, threshold resummation for the prompt photon cross 
section has been carried out by two groups || |J. Even though both results have been applied 
numerically ]10], 11], a full comparison of the two formalisms has yet to be carried out, either 
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numerically, or at an analytical level. Indeed, it should be mentioned that there are genuine 
differences between the two formalisms: [|J resums at fixed, and arbitrary, photon rapidity, 
while [[]] applies to the cross section fully integrated over rapidity. Therefore, it is interesting to 
compare the two approaches in cases where they both can be applied, which will show the close 
relation between the two. 

Providing phenomenological results for the full content of the resummation formula of H is 



our other objective in this paper. In the numerical study of the inclusive cross section in |TIJ, the 
full resummed expression of , which is at next-to-leading logarithmic (NLL) accuracy, was used, 
although rapidity acceptance was implemented indirectly, by adopting the angular dependence 
of NLO. In contrast, the analysis of [II], which made use of the formalism of ||, employed 
an expansion in the strong coupling truncated at next-to- next-to-leading order (NNLO). It is 
therefore interesting to study the role of higher-order terms, which do seem to have a fairly 
significant influence on the final result. To see how, we must develop a new practical technique 
for treating rapidity dependence in the resummed cross section, to which we will turn first in 
Sec. 2. This is an essential extension of the threshold resummation formalism, applicable to any 
resummed cross section at measured rapidity. Again, to make direct contact with the numerical 
studies of [HJ, we adopt the minimal prescription for moment inversion of fl6 |. 

In Sec. 3, we give explicit expressions for the resummation exponents at NLL accuracy, arising 
from the formalism of j8|. Section 4 presents the NLO hard-scattering functions that multiply 
the resummed exponents and collect contributions that are constant near threshold. The actual 
comparison between the two resummation formalisms of || and || is performed at an analytic 
level in Sec. 5. Section 6 discusses our treatment of the inverse Mellin and Fourier transforms 
required to obtain the physical cross section from the resummation formalism laid out in Sec. 2. 
Finally, in Sec. 7 we present numerical results for both resummation formalisms and provide 
phenomenological predictions for the prompt photon cross section resummed at fixed rapidity. 



2 The resummed prompt-photon cross section 

The cross section p^da/dprdr] for the single-inclusive production of a photon with transverse 
momentum px and pseudorapidity r\ is given in terms of the partonic cross section p\da /dpTdrj 
and the distribution functions (f) a & for initial parton types a, b as 



dpxdrj J x ™ in Jxf in dpxd-q 



where Pa-, Pb, P-y are the momenta of the initial hadrons and the photon, and 



x T e v min x a x T e v 



min = mm = = , , . 

2-x T e- r ? ' b 2x a -x T ei ' y /v v ' 

In Eq. (]]]), \i F and fj,R denote the factorization and renormalization scales, respectively. 
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Threshold resummation is organized in Mellin- iV moment space, where both the evolution of 
the parton densities and their convolution with the hard subprocess cross sections simplify. Were 
we interested simply in the cross section da/dpx, integrated over all rapidities, an appropriate 
variable for taking moments would be x\, as shown in ref. ||. If we wish to retain the rapidity 
dependence of the single-inclusive photon cross section da/dprdrj, we need to introduce a double 
transform of the cross section. It turns out to be particularly convenient to do this by defining 

^^jy^r'-^r 1 ^- (3) 



2n J-oo Jo ^ 1 dp T dr] 

In this expression, we take Mellin moments in x\ and a Fourier transform in rapidity. The latter 
makes use of the relation between the hadron and parton level center-of-mass system (c.m.s.) 
rapidities, r\ and fj, respectively: 

,-« + ih(!), (4) 

where x a ,Xb are the parton momentum fractions. Then, the Fourier phase factor exp(iMrj) 
carries through to the parton level as exp(iMrj) x % ^ [ l 2 x h lM ^ 2 j so that we will eventually arrive 
at a Fourier transform in partonic c.m.s. rapidity, along with shifts by ±zM/2 in the complex 
plane of the Mellin moments of the parton distributions. Introducing the invariant mass of the 
unobserved hadronic final state, 

S4 = s (1 — xt cosh ff) , (5) 



where 



we indeed obtain 



x a x b S , x T = XT , (6) 



a{N , M) = J 2 - •£ r (cosh,))-" f ( 1 - £i) B*%* 

V 7r TT J-oo Jo s V s J apTdn 



o,6 

r ^e^(coshr})~ 2JV ^(2iV,r)) . (7) 

271 a,b J ~°° 

The second relation serves as the definition of the hard-scattering function u) &, in Mellin- moment 
space, as a function of rapidity. Moments of the parton distributions are denoted by 

i 

dx x L ~ l (j) a (x,fj%?) . (8) 

o 

Threshold resummation for the single-particle inclusive cross section organizes corrections as 
singular as a" [ln m (s 4 /s) /s^],, m < 2n — 1, near partonic threshold, S4 — > 0, to all orders in 
perturbation theory. Let us define by N' = 2N the moment variable conjugate to (1 — s 4 /s). 
This variable was denoted as iV in Ref. ||. In Mellin- N' space, these logarithms translate into 
terms of the form a"ln m (A^'), m < 2n. The hard scattering functions u a b(N',fj) for qg and qq 
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scattering, resummed to NLL accuracy at fixed rapidity fj, were derived in ||, and are given as 



-a,b 



tiab(N',f)) = aa s (/i R ) C ab (f),a B (jj? R )) exp {E' C {N', s)} exp <j £ E^N',, s) 

S/N '^2Rer^ (fj,a s (^ 
A 4 



x exp 



(9) 



with N[ = 2Ni, N a = N(—u/s), N b = N(—t/s). We now describe each of these functions, starting 
from the left. The coefficients C a b(fj, a s ) include the underlying Born 2 — > 2 hard scatterings 
ab — ► jc, and also absorb factors in the cross section that are constant in the large- N' limit at 
fixed fj, produced by hard virtual and non-logarithmic soft higher-order contributions 0. The 
C a b are subject to 2 — >• 2 kinematics, with xt — 1/ cosh(r/), and are functions of fj and the strong 
coupling only, 

C ab (v,a s ) = f:^yC%\f l ). (10) 



Their leading terms are just the 2-^2 cross sections for ab —>■ 7c, given in our normalization 
by 



(V) 

(v) 



Tre 2 q C F 1 ft u 
C A 2 cosh 4 (77) \u + l 
Tie 2 1 / t s 
2C A 2cosh 4 (r)) V s ~ t 



(11) 



where t = (p a — p 7 ) 2 = — sx^-e *?/2, w = (p& — p 7 ) 2 = —sxt^/2. The next-order terms C„J can 



be extracted explicitly || from the complete NLO calculation of |fL7|| ; see below. 

To the right of the coefficient functions C ab in Eq. (0), each of the exponential factors is 
associated with one of the sets of on-shell quanta that characterize the hard-scattering process 
near partonic threshold. These are: quanta collinear to the incoming partons, quanta collinear 
to the final-state jet that recoils against the photon, and soft-gluon radiation. 

The resummed exponents for the initial state partons in Eq. (0) are given by 



Ei(2Ni, Mi 



,2AT;-1 



dz- 







z 



dt 



(l-z) 2 t 



A, 



a. 



{tMf) 



B i [v i ,^l,a a ((l-zyM?) 



M i dy! 



-2/ 4 7, (« s (/i' 2 )) 



dp! 



2 1 ^ lii (2N h a s {^ 



(12) 



where Mj is a scale of order ^/i. As we shall see, the exponent is actually independent of Mj to 
NLL. In addition, 



Ma s ) = Ci\ — + ^K 
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K = C / 



'67 

18 6 



n/ , 



5 ( M ? N 



7T 



/M 2N 

1 - \n(2u t ) +ln f — - 



(13) 
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Here, C q = Cp = 4/3 and C g = Ca = 3, and n/ is the number of flavors. The Ui contain the 
gauge dependence, introduced to define the factorized cross section ||, 

\ n I 

where /3f = $J2/s, and n an axial gauge vector. Furthermore, to one loop order 
7,(a.) = |C^, 7 w (^,a.) = -(lii^-|)c7 



3 „ a., , „ , /, 3\ ^ a s 

F , 

7T 

OL 

j g {a s ) = b a s , J gg (N, a s ) = - (CUlniV - nb ) — , (15) 



7T 



where 6 = (HCU - 4T R ra / )/(12vr), with Tr = 1/2, and where iV = No 1E , with 7^ the Euler 
constant. The 7« are the anomalous dimensions of the quark and gluon fields, and the 7^ the In iV 
and constant terms in the negative of the moments of the diagonal splitting functions. As can be 
seen from (|^), they organize renormalization- and factorization-scale dependence, respectively. 

The four terms in the initial-state exponent are closely linked by the relation between the 
large- iV behavior of diagonal splitting functions and anomalous dimensions, 

-fa(N,a s ) = -lnNAi(a s )+ji(a s ), (16) 

where, as above, iV = Ne lE . Each of the four terms in Ei plays a specific role in the resummation 
process. The first, double-integral, term organizes double-logarithmic behavior associated with 
soft and collinear divergences in terms of the universal behavior of splitting functions. The 
second, single-logarithmic, term matches the behavior of the splitting functions to that of the 
physical cross section, factorized as in Ref. ||. The third term summarizes noneikonal logarithms 
in virtual diagrams, while the fourth matches large transverse momentum subtractions associated 
with the MS scheme to the double-logarithmic term. Using Eq. (|nSD, we readily show that at NLL, 
Ei is independent of the Mj. In our notation, Refs. M and [nj adopt the choices M a = —uj^fs 
and M b = -t/y/s. 



From the final-state jet that recoils against the photon, one has the exponent 

'(l-z) 2 t 



z) dt 

E' c (2N,s) = I dz~— -{I -A c [a s (ts)} - lc [a s ((l - z)s) 




-B, 



i 1 rVs du! 
v c , 1, a,((l - zfs) + 2 / J^ lc (a s ^' 2 )) , (17) 



where B c is defined in Eq. (^) above. As for the initial-state jet, the A c and B c terms may be 
associated with the eikonal approximation, while the 7 C terms organize the remaining, collinear 
logarithms. Again, compared to Ref. [Q, we make the dependence on the renormalization scale 
explicit, although it is not associated with logarithms of the moment variable. 
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The final factor in Eq. (||]) resums NLL terms from coherent soft gluon emission between the 
jets. The rapidity-dependent anomalous dimension governing this soft radiation depends on the 
flavor combination in the hard scattering. For the process qq — > 7*7, it is given by 



2tt 



C F [- \n(Av q Vq) + 2 - 2vri] + C A 



In 



For the process qg — > jq, it reads 



^(99^79) 



2^ 



C h 



2 In 



-u 



-MH^9c) + 2 



ttt\ 



hi 



+ 1 - ln(2i/ fl ) + iii 
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- ln(2f s ) + 1 - 7U 



.(19) 



Note that the dependence on the gauge vector n in all these expressions, contained in the 
Eq. flLf), drops out when all contributions to Eq. (H) are combined [§, [TTj, pL8| . Indeed, by 
redefining the soft and jet functions as in Refs. we may eliminate this dependence at the 
beginning. Here, however, we retain it, again for ease of comparison to Refs. || and |JTT . 



3 NLL exponents 



Following the procedures of |], 20 1, we reexpress the three exponents in Eq. ([]) in a form with 
NLL accuracy. The respective results for each exponent are given in Appendix A. The sources 
of 17 dependence are readily visible from the expressions given there and are: the dependence of 
the initial-state exponents Ei on Aj = a s bo\iaNi (where a s = a s (fi 2 R )), the dependence of the 
final-state jet function on s = cosh 2 17, and the dependence of the soft anomalous dimensions 
p(a6^7c) Qn ^ e subprocess Mandelstam variables s, t, w, as given above. After combining all NLL 
exponents given in Appendix B, we find that for both partonic channels the 17 dependence is 
purely of the form ln(coshi7): 



[a a* 



T 1 / drje iMf > (coshfi)- 2N uj ab (2N,fj) 

J— 00 

dfj e lMf > (coshf))~ C afe (i7,a s )exp 

00 

: exp [lniV^ ) (A)+^ ) ( A, %, %) 



(20) 



InA^(A) + 9$ (a, ln(coshr)) ^(A) 

V r"R f^F / 



dfje iMfi {cosh fj) 
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where Q = 2py, and N = N+l, A = a s bo ln(iV+l). Here we have adopted the convention of || 
to evaluate the resummation exponents at the moment variable N + l, rather than at N, since 
the former naturally appears in the moments of the parton densities; see Eq. (|7|) and Eqs. (0), 
(|33| ) below. The functions g^ b \\) an d g^b \\ Q 2 1 f^w Q 2 1 \*f) (where A = a s b \nN) resum LL 
and NLL terms, respectively, that are independent of 17. Resummation of the 17 dependence is 
taken care of by the functions g^(X), which are of NLL order only. The particular simplicity 

(3) 

of the fj dependence in the exponent makes it possible to join all g ab terms with the moment 
variable in the fj integrand. 
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The functions g^ b \\) and g^{\ Q 2 /fi 2 R , Q 2 /n 2 p) have already been presented in ||, and we 
recall them here: 



gV(\) = (2C F -C A )h^(\) + C A h (1 \\/2) 
9$M = C A h^(X) + C F h( 1 \X/2), 



where 



while, 



( 2 ) (\ 91 91 

yqq I ^' ,,2 ' ,,2 

\ H'R H-F, 



1 



2tt6 A 



2A 



1 - 2A) ln(l - 2A) 



{2C F - C A ) /i (2) (A) + 2C A /i (2) (A/2) 



(21) 
(22) 

(23) 



+ In 2 ln(l - 2A) + C ™~ **» ln(l - A) - ^ A In % 



+ 



nbn 



2A + ln(l -2A) 



+ 



C A 
2nbn 



21n(l - A) - ln(l - 2A) 



In 



Q 2 



i > 



( 2) i x 91 91 
9qa[ '/4'/4, 



C A h^ 2 \\) + 2C F h { - 2 \\/2) 



(24) 



i n1 /-, OXA . 4CV 7s — 3C F C f + C a Q 2 
+ — -In21n(l-2A) + — ln(l-A) Aln — 

2tto Attoq 7ro Hp 



+ 



Op + C j. 



27rfe n 



2A + ln(l - 2A) 



27rfe n 



21n(l - A) - ln(l - 2A) 



1 

In — , 

/4 



where Q 2 = 2p T , 7^ is the Euler constant, -yp = 0.5772 . . ., and 



2nbl 



2A + ln(l-2A) + -ln 2 (l-2A) 



IE 

nb 



ln(l-2A) 



K 



4n 2 b 2 



2A + ln(l-2A) , (25) 



with bi = {17C\ — 10 C A T R rif — 6 C F T R rif) / (24ti 2 ) . The coefficient K has already been defined 
in fll3"|) . The functions g^ b \\) are new to this analysis and read 



where 



gif(X) = (2C F -C A )h®(\) + C A h®(\/2), 
4? (A) = C A h^\\) + Cph^\\/2), 



h®(\) = — W1-2A) . 

TT&o 



(26) 
(27) 



As mentioned earlier, they only contribute at NLL level. 

Finally, following Eq. (PD, we have expanded the hard-scattering functions C a b(f),a s ) in 
Eq. ( |2!1D to first order in a s . As was pointed out in || [TI|, the inclusion of the terms 
is significant, since with their help three towers of logarithms, a™ln m (iV), 2n — 2 < m < 2n, 
are brought under control in Eqs. ([]), (PU|). In this way, the factors C a b also help reduce the 
dependence of the cross section on the renormalization scale fiR- 
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4 The hard-scattering functions 



Expanding the expressions above in a s , and comparing them to the exact NLO calculation 
of [ IT | , we are able to check that indeed the logarithms of the type a s ln 2 (iV) and a s ln(iV), which 



appear at NLO, are correctly reproduced by the resummation calculation, including their rapidity 



dependence |TT| . This comparison also determines the explicit expressions for the iV-independent 
coefficients (fj). For the qq annihilation process, we find: 

Q 2 

^qq' { r D S \4 1II[V{1 — V)\yi ■+■ 111 ['1^1 — "ill — \n 111 - 

20 

y 



2p 2 F v(l - v) 



-41nvln(l -v) + ^n 2 + p\ q + 6p qq - - 19 



C A 12(1 + Pqq ) \n[v(l - v)] + 24 \nv ln(l - v) + tt 2 + 3p 2 g - + 18p qq - 5 



+8vr6 29 + 3p qq - - 12 In 



Q 2 



2p 2 R Jv(l - v) 



7re 2 C F 

H - v 

C A 



l-v){ - (C A -2C F )(l + 2v) \n 2 v-2C A \nv + 2C F (5-6v) In-, 



+ (v < — > 1 — v) 



(2f 



where 



! t e 5 ? 

s 2 cosh 17 
-3 + 4 7i; + 21n[4t;(l - v)] 



(29) 



For QCD Compton scattering, one has: 



\ C<S (V) \C F [4 In 2 + 2(ln v - 3) In v - p%> In 



Q 2 



2p 2 F (l-v) 
1 / (p\\ 2 3 fp * 11 2 29\ 
+ 8 ("« ) + 2 P "» + ¥ _ T j 



C ( In 2 — + In 



2/i^ 2 



7re z 

+^ 2 (1-^)|(^-2C F ) 



1-v 



(3 - 2v) In 2 

V v 



+ tt 2 + (1 - 2v)\n 2 v 



+2(1 -u) In 



1-v 



+ 6C F \n(l-v) 



(30) 



where 



P (F) 

P (A) 
rqg 



-3 + 4( 7£; + ln[2' 
4( 7i? + ln[2(l -*;)]) . 



(31) 
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The values of the ratios (rj) / (rj) at rj = (that is, at v — 1/2) agree with the corresponding 
expressions in ||. Despite the somewhat complicated rapidity dependence of the C^ h \ it is 
possible to evaluate the integral in the last line of Eq. (^) in closed form. These results are given 
in Appendix B and complete the derivation of the resummed photon cross section in Mellm-TV 
and Fourier space. We note that the explicit 7£-terms of the coefficients, in Eqs. ( |2"§D and (|3~T|), 

fo\ 

as well as in g ab in Eqs. (|23J) and (^), can be eliminated [jnj by redefining A = a s bo (lniVe 7J5 ) = 
a s bo\nN. This convention was adopted in Ref. f!M. 



5 The rapidity-integrated resummed cross section 



We now compare to the calculation presented in ||, which aimed at threshold resummation of 
the rj- integrated cross section. According to Eq. (Q), integration over r/ is equivalent to setting 
M = in our expressions and multiplying by \^2tt. We then find from (|3]), (|7]), ( p0|) for the 
resummed cross section: 



rl/ cosh 2 rj 

di] / dx 
Jo 



2 \N-i p^da 



(4) 



a a. 



E 

0,6 



^ +1 0f +1 ex P 



dpTdrj 



(32) 



dfj (cosh 17) 2N+9 ' 



(3) 



(A) 



at. 



1(1), 



7T 



to be compared to the result of p[, which can be written as 

N-i pj,da 



oo rl/ cosh 2 rj 

drj \ 

-oo Jo 



dcc^p ( oc^p 



0,6 



" +1 ^ +1 exp 



dp^d?7 
lniV^W+^W 



(33) 



«.ci?(o) 



1 + 



7T 



C2>(0) 



dry (coshr}) 2JV C^r?) 



The two results become identical if the resummed exponent and the first-order correction in 
C a b in Eq. ( [32]) are evaluated at 77 = 0, instead of retaining their full rj dependence. We note 
from (0), however, that for large N the integrand is indeed strongly dominated by the region 

rj ~ 0. Therefore, setting C^(rj) / 'C^ h (fj) = (0)/C^° (0) and neglecting the term proportional 

(3) 

to g ab is expected to be a good approximation, and we anticipate that the two resummation 
formalisms should yield numerically very similar results for the //-integrated resummed cross 
section. On the other hand, if one wants to obtain fully consistent resummed results at fixed 
rapidity, one needs to make use of the full expression in Eq. (|20| ) resulting from the formalism 
of i. 
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6 Inverse transforms 



The final step is to insert the result in Eq. ( |2T)[ ) back into (|7|), and to take the Mellin and Fourier 
inverse transforms of <r(N, M), in order to arrive at the physical hadronic cross section: 



p^dcr 
dprdrj 



1= r dMe- iM ^ I 
2n J-oo 2m Jc 



C+ioo 



V2vr 
aa s 
~2t7 

x exp 



dN (xl 



-N 



oo 
oo 



dMe -iM V 
oo 2m JC-ioo 



C—ioo 
C+ioo 



a{N, M) 



dN (xl 



-n i N + 1 + 1 ¥- i N+1 - l -f- 



a.b 



\nNg^(X)+g ( j ! 
dr)e iMf > (cosh fj) 



\ 91-9^ 

' ;; 2 ' ii 2 

\ H'R PF, 



-2N+g$>(\) 
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(34) 



where for the second equality we have made use of Eqs. (0) and (p0|), and where the last integral 
is given explicitly in Appendix B, as discussed above. Care has to be taken when choosing the 
contour in complex N space. First of all, the functions ^(A) have cut singularities starting at 
A = J; that is, at iV = N L = exp(l/2b a s (fi R )). These singularities result from the sensitivity 
of the original resummed expression (P) to the Landau pole [20j, and signal the onset of non- 



perturbative phenomena very close to threshold. With the "minimal prescription" [[K| |16j for the 
exponents given in Eqs. ([0])-(26D above, we choose the constant C in (0) so that all singularities 
in the integrand are to the left of the integration contour, except for the Landau singularity at 
N = Nl, which lies to the far right. The contour is then deformed [10, 16, [H[ into the half-plane 
with negative real part, which improves convergence while retaining the perturbative expansion. 
In this deformation, we need to avoid the moment-space singularities of the parton densities, 
which are displaced parallel to the imaginary axis by ±M/2, as seen from Eq. (^). Thus, the 
intersection at C of the contour with the real axis has to lie far enough to the right that the 
contour does not pass through or below the singularities of the parton densities. At the same 
time, one can make the bending angle k close to 90° to prevent this from happening. Our 
parameterization of the contour is given as 



N 



C + ze 



7r — arctan 



(0 < z < 00) , 
' C-l + M/2 \ 
, C-l J 



(35) 



The situation is depicted in Fig. [I]. 



7 Numerical results 



For our numerical calculations, we use the NLO parton densities of [22]. The evolution code 
of [E2] is set up in Mellin- N moment space and is therefore ideally suited for implementing the 
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Figure 1: Contour in Mellin-iV space for inverting the resummed cross section. The asterisks 
denote the rightmost poles of the parton densities which have acquired an imaginary part through 
the Fourier transformation in rapidity. Nl is the position of the Landau pole; see text. 



resummation formulas in the way we have outlined above. In order to avoid double-counting 
of higher-order corrections, we subtract from the resummed expression (|9|) its perturbative ex- 
pansion to 0(a s ), and subsequently add back in the full "exact" NLO prompt photon cross 
section ITO] . For the latter, we use the calculation and the numerical code of [Ol and also include 



an NLO fragmentation component, computed using the results of [p3fl , along with the NLO pho- 
ton fragmentation functions of [24]. Note that complete consistency would demand to also take 
into account the effects of threshold resummation on the fragmentation component [10, [14]]. To 
mimic the expected effects, we multiply the NLO fragmentation component by the same factor 
by which the non-fragmentation part is changed through threshold resummation. Fortunately, 
the fragmentation contribution to the cross section is generally at the level of only 15 — 20% 
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anyway || |10|, so the uncertainty residing here is not really crucial. 



As in [pi 11|, we will focus our calculations on the fixed-target region, where the observed 



disagreement of data and NLO theory calculations is largest, and which also is known [ID, |TT|| to 
be generally more susceptible to threshold resummation effects than the collider regime, where the 
relevant values of xt are much smaller. Figure ^| compares our results for the dependence of the 
resummed prompt photon cross section with those of |]IOfl . We use the same parton distribution 
functions [22] and rescaled fragmentation contribution for both curves^. We show E^cPa/dp^, 
calculated for p-N collisions at -Ebeam = 530 GeV, corresponding to the Fermilab E706 fixed- 
target experiment 0, whose data are also displayed. We have chosen the scales Hf = = Vt 
and integrated over the rapidity range \r}\ < 0.75 accessed in the E706 experiment ||. For better 
visibility, we have normalized all cross sections by the NLO one. One notices first of all that 
threshold resummation induces large effects at high p^, where it leads to a sizeable enhancement 
of the theory prediction [II], |TTJ , but is unimportant at the low-pr end. Figure ^ also shows that 



there are only small differences between the results given by the two resummation formalisms. 
These have two origins: first, they are related to the neglect of the rapidity dependence in 
the resummation exponent and in the higher-order terms of the iV-independent coefficients C a b 
in [[J. As mentioned above, these differences would even be there, had we integrated over all 
kinematically allowed 77, that is, over \rj\ < In ( 1 + Jl — x\ ) fx? ■ Second, since rapidity 



1 + ^ji-xiyxT 

dependence in resummation was not available in |H|, the means to compare to experimental 
data covering only a limited range of rapidity was to apply the following acceptance correction 
to the resummed cross section: 

„<~>(M < 0.75) S ,(»)(all ,) "'""'^lall") 75 ' ' (36 > 

That is, it was assumed that rescaling the //-integrated resummed cross section by an appropriate 
ratio of NLO cross sections correctly reproduces its rapidity dependence. This assumption of 
similarity of the resummed and NLO ?7-shapes is actually a fairly good approximation, as we 
will see below. We conclude from Fig. |2| that retaining the rapidity dependence in resummation 
induces very mild overall corrections for a cross section that is integrated over a substantial region 
including rj = 0. As is evident from Figure threshold resummation effects, albeit sizable in 
certain kinematic regions, are not sufficient to obtain a satisfactory description of the E706 data. 
This conclusion was first drawn in Refs. [TO] , [TT |. 



We now turn to the rapidity dependence of the threshold-resummed cross section, where using 
the full formalism of || presented above is mandatory. The UA6 collaboration has presented 
data for the prompt photon cross section as a function of rj, taken in pp and pp fixed-target 
scattering at v^S* = 24.3 GeV M. The transverse momentum of the photon was averaged over 
4.1 < p T < 7.7 GeV. In Fig. |3| we compare our calculations to these data. In both plots, 

lr The precise treatment of fragmentation and the choice of parton densities is different than in [lQ[| . These 
choices have little numerical effect, however. 
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Figure 2: Threshold-resummed prompt photon cross sections, normalized to NLO, based on the 
formalisms of |J (dashed) and || (solid). The parton densities were taken from [ 2"2"| . The E706 
prompt photon data H are also shown in the same normalization. 



the dashed lines display the fixed-order NLO cross sections for three different choices of the 
renormalization and factorizations scales, fip — — 2Pt,Pt,Pt/2, from lower to upper. As can 
be seen, the variation of the NLO predictions induced by scale changes is very large, indicating 
the importance of yet higher-order effects. The solid lines in Fig. |^ show the results after adding 
in the effects of NLL threshold resummation. The most striking feature is a strong reduction 
of scale dependence. This effect was also found for the 77-integrated cross section in our 
results demonstrate that it occurs over the whole 77 range. A closer analysis reveals that the 
main reduction in scale dependence is associated with the factorization scale. The reasons for 
this were discussed in [13|, |25fl . Also, the scale dependence decreases even further if we push our 
calculation to larger xt- Both these observations are corroborated by Fig. f|, where we return to 
the conditions of the E706 experiment ||, looking at the rapidity dependence of the cross section 
at a fixed, high, photon transverse momentum px = 10 GeV. The dashed lines again show the 
NLO cross section for three choices of scales, \ip = fiR = 2pt,Pt,Pt/2, and the solid lines display 
their resummed counterparts. Clearly, the reduction in scale dependence after resummation is 
dramatic here. The dotted lines in Fig. ^ refer to the NLO cross section for fixed renormalization 
scale hr = pt, but varying the factorization scale in the range used previously. The scale 
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dependence of the NLO cross section thus mainly is due to that on //p, and the improvement in 
fip dependence after resummation, shown by the dot-dashed lines, is striking. 

We note that resummation of the ^-dependence of the cross section was also considered in [IT|. 
There the resummed expression was expanded to O(o%), thus testing the effects of the NNLO 
terms generated by resummation. The numerical results of |TTJ indeed also show some reduction 
in scale dependence which, however, is not as substantial as the one we find. This implies 
that resummation effects beyond NNLO are still important, and that the full resummed series 
should be taken into account. Our examples in Figs. |3|, |] also show that threshold resummation 
mainly affects the normalization of the cross section, and less so its shape in rj for moderate rj. 
Nevertheless, as we now show, large rapidity can, like high px, enhance resummation effects. 




Figure 3: Rapidity dependence of the prompt photon cross section in pp and pp collisions at 
y/S = 24.3 GeV. Dashed lines are NLO, solid lines denote the cross section resummed to NLL 
accuracy. For each case, the results have been calculated for three choices of scales, fip = = 
2Pt,Pt,Pt /2, from lower to upper. The data are from UA6 HQ. 
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Figure 4: Rapidity dependence of the prompt photon cross section in pN collisions at = 31.5 
GeV. Dashed and dotted lines are NLO, solid and dot-dashed lines denote the cross section 
resummed to NLL accuracy. For each case, the results have been calculated for five choices of 
scales, fip — [Mr — 2p T ,p T ,p T /2 (dashed and solid), and fip = 2p T ,p T /2 at fixed = pp (dotted 
and dot-dashed), from lower to upper. 



We have mentioned earlier that threshold resummation for the prompt photon cross section 
is typically expected to be relevant in the fixed-target regime, but less so for collider energies, 
where a given xp implies a much larger value of pt and hence a much smaller cross section. As 
a result, the highest-pr prompt photons seen at the Tevatron are at around pp ~ 100 GeV 0, 
where iy pa 0.1 which is too small for threshold effects to be really important. However, we 
can ask what the effects of threshold resummation would be, were we able to go out in pp 



to, say, pt = 500 GeV. Keeping in mind that the Tevatron jet data [|(J do reach out to pp's 
of that size, this question becomes of some relevance, provided threshold resummation effects 
for the prompt photon cross section have any resemblance to those for the jet cross section. 
Fig. |^ shows a calculation of threshold resummation effects on the cross section for very-high-pr 
prompt photons, produced in pp collisions at \fS = 1800 GeV. The results are normalized to the 
corresponding NLO cross section. The solid line shows the cross section averaged over |?7| < 0.9, 
used for the actual CDF prompt photon data 0, whereas the dashed and dotted lines present the 
result at a fixed near-central (t] = 0.7) and forward (r] = 1.4) rapidity, respectively. It is evident 
that at such high pt threshold resummation effects are sizable even at collider energy, and that, 
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at a given p?, they tend to become more important with increasing \rj\, that is, toward the limit 
of phase space where large momentum fractions in the parton densities are probed. At low rj, 
effects remain modest even at large p?, as recently observed directly in an NNLO expansion [27| 
of the resummed jet cross section. Our finding could be indicative of significant effects also for jet 
production at the Tevatron, in particular in the off-central region. This result is also suggestive 
of recent data PJ which hints at an excess of high-p r jets that increases with 77. This issue will 
bear further investigation. 




Figure 5: Threshold-resummed prompt photon cross sections, normalized to NLO, for pp colli- 
sions at VS = 1800 GeV. 



Acknowledgments 

We are grateful to S. Catani and N. Kidonakis for very valuable discussions. This work was 
supported in part by the National Science Foundation, grant PHY9722101. G.S. acknowledges 
the hospitality of Brookhaven National Laboratory. W.V. thanks RIKEN, Brookhaven National 
Laboratory and the U.S. Department of Energy (contract number DE-AC02-98CH10886) for 
providing the facilities essential for the completion of this work. 



17 



Appendix A 



In this appendix we give explicit expansions of the exponents in Eq. (^) to NLL accuracy. These 
are analogous to the expressions given in Appendix A of Ref . 0] . We expand each of the functions 

Qi{a s ) = Ai(a s ), Bi(. . . ,a s ), ji(a 8 ), ju(a 8 ), ReT^" 1 ^ (fj,a s ) as a power series in a s /7r: 



fc=i v ^ 



(37) 



We then find 

Ei(2Ni,Mi 

s) = 



-H[ (A*) + 7i • ; I Aj, — |-, — |- 

1 \ ^R ^F / 



M? MT 



B?> 1^,^1^(1-2^) 



M 



2nb c 



^ 0) (A)+^(A, 



a s (jj, 



(i> 



Sf } (^,l)ln(l-2A) 



2nb ( 



I'R 



— 2ReT 



(ab— >7c) 



r2- 



,(a&-> 7 c)(l) 



ln(l - 2A) 



(3? 



where A = a s (fi R )bo In N, Aj = a s (fi R )b In iVj, and 



A 



^ (0) (A) 

' 2 ' 2 / 
^R J 



^ (0) (A) 
-^ (1) (A) 



.4 



(i) 



27T&2 

2vr6[j 



+ 



2A + (1 -2A)ln(l -2A) 

- ln 2 (l - 2A) + 2A + ln(l - 2A) 



(39) 



Af ) ( 7E + ln2) 



nb 



ln(l - 2A) 



A 



(2) 



2nbn 



7ib 



2Hf\\/2)-rif\\) 



+ AP In 



(0), 



,/4 



2r^\\/2)-^\\) 



A?\ lE + \n2)-^> ln(l-A) 



2A + ln(l - 2A) 



Ml 



A- 



(i) 



nb n 



-A In 



,/4 



(40) 



To NLL accuracy, we further expand 



ln(l - 2A,) » ln(l - 2A) - H^/N) 



(41) 



where we recall that N a /N = —u/s, N^/N = —t/s. From the appropriate sums of these ex- 
pressions, the ^'s of Eqs. (|2"I]) -(p7D are identified through the linear dependence on IniV and 
In (cosh fj). 
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Appendix B 



In this appendix, we exhibit the ^-integrals of Eq. ( PQ ) for the coefficients Cj^'^^fj) of Eqs. (0), 
fl3U|). After substitution x = exp(2r/), one is led to integrals of the type 



POO 

P cd (a,b) = 2- 2 - b dxx a (l + x) b \n c x\n d (l + x) (c, d G {0, 1, 2}) , 

JO 

which can be expressed in terms of the T function and its derivatives: 

i /i b\ r(o + i)r(-a-6-i) 
Po ' o(G ' 6) = 2 5 [r "2 J f^-pj ' 

Pi,o(a,6) = P 0) oM) [#» + 1) -^(-a-6- 1)] , 

P 0) i(a,6) = p 0j0 ( a , 6) [^(-6) - -0(-o - 6 - 1)] , 

P 2 ,o(a,6) = P 0i0 (a,6) [^(a + l)+^(-a-6-l)] + [P li0 (a,6)] 2 /Po,o(«,fe) 

P 0)2 (a, 6) = P 0i0 (a, 6) [^(-a - 6 - 1) - ^(-6)] + [P Q) i(a, 6)] 2 /P 0)0 (a, 6) , 

Pi,i(a,6) = P , (a,6) ip'(-a-b- 1) + Pi, (a, 6)P ,i(a, 6)/P 0i0 (a, 6) , 



(42) 



(43) 



where ip{z) = dlnT(z)/dz and ip'(z) = dip(z)/dz. The numerical treatment of these functions 
for complex arguments was discussed in |[28|| . For convenience, we also abbreviate 



Pi j (a, 6) = Pj(a, 6) + P^a + 2, 6) , 

S M (a, 6) = 2P M (a, 6) + 2P ij (a + 1, 6) + Pj(a + 2, 6) . 



(44) 



We then have for qq scattering: 



^e^(coshr})- 2Ar+ ^ )(A) 



71 



C F 7T 2 



e 2 q {R 0fi (a, b) 1 + 



o:. s 



7T 



(2C F - CU/2) ( 7 I + 2C(2)) + 7Tb ( 7E + 2 In 2) 



+ (6C F - C a ) 7e In 2 - ~C F In 2 + (4C F - CU/2) In 2 2 
+K/2 — K q — C F {2 lE + 2 In 2 - 3/2) In \ - tt6 In \ 



{Q 1 "1 

- -C F Po,i(«, 6) - T(P 2 ,o(a, b) + 3P , 2 ) (3C A - 10C F ) + -Pi,i(a, 6) (5C A - 1GC F ) 

+ (p li0 (a, 6) - 2P 0jl (a, b)) [~\c F + (4^f - C A ) lE + tt6 - C F In ^ + (5C F - C A ) In 2^ 



+^C F P 1>0 (a, 6) + ~(C A - 2C F ) ( - 2P li0 (a + 1,6)+ 4P 0)1 (a + 1,6) + P 2 , (a, 6) 
-2P 2 , (a + 1, b) - 4P , 2 (a +1,6)- 2P 1 , 1 (a, 6) + 4P ljl (a + 1, 6) 



(45) 
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where 



K, 



\ ~ C(2)) C F , 

a = N-g®(\)/2 + iM/2, 
b = -2N + g$(\)-4. 

The coefficient K has already been defined in Eq. (|13|) . For qg scattering, we find: 



(46) 



— oo 
71 



df)e iMfl (cosh fj)' 



(47) 



Of, 



e { S 0fi (a,b) 1 + — 



7T 



(CU + C F /2) T ! - + (CU + 2C F )C(2) + (2CU + 3C F /2) In 2 2 



+ (3C A + 2CV) 7f In 2 + vr6 In 4 + 7^ f 7s + In - (CU + C F ) ( 7f + In 2) In ) 
+^|s , i,o(a, 6) (-^C F + (2a + C f ) 7f + 2(Oi + C F ) In 2 - C F In ^ 
-5-0,1 ^2(2CU + C F ) lE + (5CU + 3C F ) In 2 - (Oi + C F ) In ^ + ±(CU + C F )P 0|1 (a, 6) 
-~(CU - 8C F )S 2}0 (a, b) + i(HCU + 10CV)5 , 2 (a, 6) - ~(3CU + 8C F )5 1>1 (a, 6) 



+~(C A - 2C F ) ((P ,o(a, 6) + 2P ,o(a + 1, 6))tt 2 - 4P lj0 (a, 6) - 4P 1)0 (a + 1, 6) + 2P 0)1 (a + 1, b) 
+4P 2 , (a, b) + 4P 2i0 (a + 1,6) + 3P , 2 (a, b) + 2P , 2 (a +1,6)- 6Pi,i(a, 6) - 4P ljl (a + 1, 6)) 1 1 , 



where in this case 



a = N-gf^X)/2 + l + iM/2 
b = -2N + g®(\)-5. 



(4£ 



It is evident from Eq. ( |43D that when performing the 17 integrations one inevitably introduces non- 
leading contributions in the moment variable N. This is already true when keeping only the Born 
approximations C^(fj) for the hard scattering, and is not in conflict with the approximations 
made in threshold resummation, where the resummed exponents are always convoluted with 
smooth functions, like the parton densities, that abundantly contain terms subleading in N. We 
note that for large N the one-loop parts of Eqs. fl45|) , (|47D converge to the respective coefficients 
given in Eqs. (58), (59) of §| for the case of the ^-integrated cross section, as they must. 
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